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Abstract 

OO 



Transport and mixing properties of surface currents can be detected from altimet- 
ric data by both Eulerian and Lagrangian diagnostics. In contrast with Eulerian 
^ ■ diagnostics, Lagrangian tools like the local Lyapunov exponents have the advan- 

tage of exploiting both spatial and temporal variability of the velocity field and 
are in principle able to unveil subgrid filaments generated by chaotic stirring. How- 
ever, one may wonder whether this theoretical advantage is of practical interest 
in real-data, mesoscale and submesoscale analysis, because of the uncertainties and 
resolution of altimetric products, and the non-passive nature of biogeo chemical trac- 
ers. Here we compare the ability of standard Eulerian diagnostics and the finite-size 
Lyapunov exponent in detecting instantaneaous and climatological transport and 
mixing properties. By comparing with sea-surface temperature patterns, we find 
that the two diagnostics provide similar results for slowly evolving eddies like the 
first Alboran gyre. However, the Lyapunov exponent is also able to predict the 
(sub-)mesoscale filamentary process occuring along the Algerian current and above 
the Balearic Abyssal Plain. Such filaments are also observed, with some mismatch, 
in sea-surface temperature patterns. Climatologies of Lyapunov exponents do not 
show any compact relation with other Eulerian diagnostics, unveiling a different 
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structure even at the basin scale. We conclude that filamentation dynamics can 
be detected by reprocessing available altimetric data with Lagrangian tools, giving 
insight into (sub-)mesoscale stirring processes relevant to tracer observations and 
complementing traditional Eulerian diagnostics. 

Key words: Submesoscale, filaments, altimetry, Mediterranean circulation, 
Lagrangian dynamics 



i 1 Introduction 



Satellite high-resolution daily images of tracers like sea-surface temperature and chloro- 
phyll show a large mesoscale and sub-mesoscale heterogeneity and patchiness, typically 
in the form of filaments. The full range of tracer variability observed in a climatologi- 
cal mean at the basin sc ale, can occur in da ily snapshots over spatial scales of just 10- 
100 km (see for instance iLehahn et al.l ( 20071 ) for the case of chlorophyll in the NE At- 
lantic). This large hetcrcogenity occurring on relatively short distances is able to induce 
very strong tracer gradients and to impact on important aspects of the ocean dynamics like 
lateral transport, upwclling/downwclling, and mixing. Mesoscale and sub-mesoscale variabil- 
ity is an important component of plankton dynamic s (lAbraham . 1998; Lopez et al.1 . 12001 ; 
Abraham etall 120001 iBovd et all boOOt iToner et all l2003t iMarthl l2003h . larval transport 



( Bradbury and SnelgroveL 120011 ). as well as the dispersion of contaminants and oil spills. In 

13 polar regions, (sub-)mesoscale tracer variability has been recently recognised to impact also 

14 on the thermohaline properties of the mixed layer, affecting the low frequency ocean cir - 

15 culation through the formation of waters below the mixed layer ( Sallee et all 120061 12007 ). 
i6 Mesoscale and sub-mesoscale gradients of heat and salinity can induce cells of very strong 
17 vertical velocities (several tens of meters per day) from subsurface down to and below 
is the mixed layer depth. These vertical velocities may fertilise the photic layer, creating lo- 

19 cal plankton blooms an d affecting the sequestration of organic matter from the surface 

20 ( Lew et all l200l|^ 20051). This m echanism can affect the properties of deep waters when 

21 they are formed (jPaci et all 120051 ) . Tracer gradient intensification is also a precondition to 

22 mixing. Horizontally, the elongation of water masses in thin structu res intensifies local gradi- 

23 cnts, greatly enhancing background diffusion (jLapevre et al.l . 120061 ). Vertically, the trapping 

24 of inertial waves inside filame nts creates mixing hotspots that extends to the bottom and 

25 below the mixed layer depth ( Young et all 1 19821 ). 



Tracer mes oscale and sub-mesosca l e patc h es are often found as s ociated to surface mesoscal e 
eddies (e.g. iRobinson A. R. et all Jl993l);lMcGillicuddv et ail (ll998[ );l Martin et al.l (J2002I); 
Abraham and Bowenn 20021 ) iMorrow et al.l(|2004l ); lLegal et al.l ( 2006 )h Tracers are stretched 



into filaments by the shear-dominated regions in between mesoscale eddies while the recir- 
culating regions inside eddies' cores can trap and transport tracer anomalies for timescales 
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31 comparable to the eddy lifetime. As a well-known tool to reveal geostrophic velocities asso- 

32 ciated to eddies, altimetry data when opportvmly analyzed, provide an important source of 

33 information for the mesoscalc and sub-mcsoscale structuration of tracer distributions. 
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Traditionally, transport information has been inferred from altimetry with Eulerian diag- 
nostics like the eddy kinetic energy, the Okubo- Weiss criterion that separates the eddy's core 
from the periphery and others, that will be discussed in detail shortly. These diagnostics are 
based on the analysis of instantaneous snapshots of tracer and velocity fields. For instance, 
the Okubo- Weis s (OW) parameter has been applied to drifters in ord er to characterize 
eddy properties (|Stocker and Imbergerl . 120031 : iTestor and Gascardl. 2005), to in situ densit y 



2004) 



profiles to identify regions with different mixin g properties fll scrn-Fontan et et all .. 
and t o altimetric maps with similar objectives ( Isern-Fontanet et all 120041 ; [Waugh et al 
2006) . The OW parameter has also wid ely been used for the identification and charac- 



terization of ocean eddies from altimetry ( Isern-Fontanet et all 120031 ; iMorrow et all 12004 ; 
IChaigneau and Pizarrol l2005t llsern-Fontanet et all 120061) . More recently, it has been sug- 
gested that Lagrangian diagnostics are more appropriate to link turbulence properties to 
tracer dynamics. The most used Lagrangian diagnostics have been the calculation of local 
Lyapunov exponents, that measure the relative dispersion of advected particles. In contrast 
to Eulerian tools, Lagrangian diagnostics do not analyze instantaneous snapshots of the 
velocity field, but measure transport properties along particle trajectories, therefore recon- 
structing the fine structure of transport dynamics that a fluid parcel has experienced, like 
subgrid filament formation. These patterns depend on the advection history along a tra- 
jectory, that may span a large spatiotemporal domain of the velocity field, and therefore 
cannot be captured by Eulerian diagnostics, that instead measure local, properties. Regard- 
ing Lagr angian tools like the finite-si ze or finite-time Lyapunov exponents (resp. FSLEs and 
FTLEs) lAbraham and Bowenl ( 2002 ) compute d FTLEs from time-dependent velocity fields 
and compared to istantaneous SST patterns, lLaCasce and Ohlmannl (120031) use d FSLEs 



to study the statistics of a large data set of drifting buoys; I Waugh et alj ( 20061 ) applied 
FTLEs to alt i metry and compared the results with OW and kinetic energy estimations; 



Lehahn et al.l (|20071 ) computed FSLEs from the geostrophic field in order to ex tract trans- 



port b arriers and compared with satellite derived plankton and SST patches; Ros si et al 



(|2008f ) compared mixing and productivity by using the FSLE. Besides these studies, the 
FSLE has been also applied to numerical simulations of t he Mediterranean sea with different 
objectives: the predictability of Lagrangian trajecto ries (jludicone et all 2002) and the char- 



acter ization of flow patterns (jd'Ovidio et al 
120051) . 



20041 ) and dispersion (jGarcia-Olivares et al. 
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However, for the case of real velocity data and given the chaotic nature of mesoscale advec- 
tion (in the sense of positive Lyapunov exponents) one may expect Lagrangian diagnostics 
to be more affected than Eulerian tools by any space and time indetcrmination of the ve- 
locity field. When dealing with real altimetric data, therefore one may ask whether the 
theoretical advantage of Lagrangian techniques of detecting small scale transport structure 
is not countered in practice by the intrinsic error and cutoff of altimetric data, so that this 
information is either unreliable or already contained in Eulerian quantities. This question 
is particularly relevant given the larger computational costs (a few orders of magnitude) 
of Lagrangia n tools in respect to Eulerian ones. In fact, for the case of the Tasman sea 
Waugh et al.1 ( 20061 ) found a striking resemblance between maps of Lyapunov exponents 



and the eddy kinetic energy, as well as a compact relationship between the two quantities. 
In their concluding remarks. I Waugh et al.l (J2006J) argued that this observation, if valid in 
other oceanic basins as well, would rise the possibility of using the EKE for estimating 
stirring rates, without the need of explicit, longer Lagrangian calculations. However, one 
of the main advantages of the Lagrangian techniques is the possibility of reconstructing 
tracer patterns that are below the resolution of the velocity field and that arise due to 
several iterations of stretching and folding during the tracer (chaotic) advection. Indeed, in 



83 IWaugh et al.l (|2006T ) the correlation between the Lyapunov exponents have been found on a 



B4 grid at relatively low resolution (0.5 deg.) and for a advection times of 14 days (two altimet- 

85 ric images) , i.e. without the spatial and temporal information needed for the dcvelopcmcnt 

86 of filaments. Docs this result hold also for filament-resolving Lagrangian calculations? 

87 In this paper we address this question, exploring the possibility of obtaining reliable in- 

88 formation with Lyapunov exponents that is not already contained in Eulcrian diagnostics. 

89 focusing on filament dynamics. 

90 This will be done both on the analysis of individual patterns (validated by high resolution 

91 SST images) and on a climatological basis. We find that at (sub-)mesoscale resolution the 

92 information provided by Eulerian diagnostics and Lyapunov exponents coincindes only for 

93 very stationary eddies, while providing two distinct and complementary pictures of the 

94 circulation in all the other cases: the Eulerian analysis provides the eddies that populate 

95 the mesoscale, while the Lagrangian analysis yields the tracer filaments generated by the 

96 spatiotemporal variability of these eddies. Our results also show a surprising reliability of 

97 altimetric data at the scale of their nominal resolution, when reprocessed with Lagrangian 

98 tools. 

99 We will compare Eulerian diagnostics and Lyapunov exponents in the Algerian basin because 
ioo of its rich and variable mesoscale activity, that contains jet-dominated regions (the path of 

101 the Algerian current) as well as eddies with very different characteristics (from the quasi- 

102 stationary Alboran gyres, to the slowly propagating eddies of the Balearic Abyssal plain, 

103 and the fast-evolving Algerian eddies). A main characteristic of the surface circulation in 

104 the Mediterranean sea is the propagation of fresh waters incoming from the Atlantic ocean. 

105 At the entrance of the Mediterranean, these waters flow from west to east and form patterns 

106 such as the Alboran eddies east of the Gibraltar strait or the Algerian current along the 

107 Algerian coast. The instabilities of this current generates, a few times per year, coastal 

108 eddies that propagate downstream, usually until the entrance of the Sardinia chann el. There, 

109 they can detach from the coas t and propagate as o p en se a eddies (JMillotl . [l999j) following 



relatively well defined paths 1 Isern-Fontanet et al.l . [2006). These eddies, called Algerian 



eddies, have variable diameters of about 50 — 200 km, vertical extents from hundreds to 



112 thousands of mete r s, and lifetimes of several months, up to nearly 3 years (jMillot et al 

113 



19971 ; iPuillat et al.l . 120021 ) . Their presence has a large impact on the redistribution of tracers 
in the Algerian basin, which is characterized by the nort hward spreading of tracers that are 
initia lly transported eastwards by the Algerian current ( OvchinikovL Il966c iBrasseur et al 



ii6 Il996l ). As we shall see, our analysis of stirring will unveil the crucial role on tracer patterns 

ii7 of the time variability of the mesoscale activity in the Algerian basin and the relevance of 

us topography in constraining the dynamics of coherent structures. 

ii9 The paper is organized as follows. After describing the data sets and the techniques, we de- 

120 velop an Eulerian and Lagrangian analysis for specific days. We focus on the representation 

121 of eddies and the detection of transport barriers and we get two complementary pictures: 

122 a regular and smooth OW-based description and a more complex, lobular representation 

123 from the FSLE. We then compare such different structures to tracer distributions. In order 

124 to filter out the active dynamics of a real tracer as well as the unresolved components and 

125 indctcrminacies of altimetry data, we first consider the filaments of a synthetic tracer, that 

126 we advect numerically with the altimetric data. As a second step, we take sea-surface tem- 

127 perature (SST) satellite images. The indications that are found on the individual days are 

128 then generalized in a climatological comparison, where we compare the spatial variability 

129 of FSLE, OW parameter, strain rate, and eddy kinetic energy temporally averaged over 

130 the period 1994-2004. For both the Eulerian and the Lagrangian analysis, we also propose 

131 to describe the consequences of filament formation with a climatology giving the spatial 

132 density of transport barriers. 



133 2 Methods 
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A traditional approach to the characterization of the stirring and mixing in the presence of 
mesoscale eddies consists in separating the stagnation region at the eddy's core from the 
eddy's periphery where tracer filamentation occurrs. This is done by measuring the relative 
dominance of vorticity and deformation. One of the most used parameters f or me asuring this 
relative dominance is the Okubo- Weiss parameter ( Okubol . Il970l : IWeissl. 1991) which has 



139 been already used by some of us to study properties of Algerian eddies (Jlsern-Fontanet et al 



20041 ). The OW para meter is a particul a r case of the more general vortex- identification 



criterion proposed bv ljeong and Hussainl ( 1995 ) 



142 The Okubo- Weiss parameter W is defined as: 

143 W = Sl+Sl- CJ 2 



(1) 



where s n , s s and lo are the normal and the shear components of strain, and the relative 
vorticity of the flow defined, respectively, by 



du dv 
dx dy' 



dv du 
dx dy' 



dv du 
dx dy 



147 In the formula above, x and y are orthogonal spatial coordinates and u and v are the 

148 component of the velocity respectively for the x and y directions. 
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The parameter W, allows to separate a two-dimensional flow into different regions: a 
vorticity-dominated region (W < —Wo), a strain-dominated region (W > Wo) and a 
background field with small positive and negative values of W (\W\ < Wo). Here Wo = 
0.2ctw-, <jw being the stan dard deviation of the W values in the whole domain, in our 
case the Mediterranean sea (JBracco et all l2000t IPasquero et all l2001t llsern-Fontanet et al. , 
2006; El hmaidi et all Il993l ). The core edge can then be identified as the closed lines with 
W = 0. This separation of the field in terms of the sign of W has been proved to be a 



156 robust criterion for extract ing eddy cores from complex fluid flows (jJeong and Hussain 



1995c IPasquero et all l200lh . In steady flows, the boundary of the core constitutes a bar- 



rier to the exchange of particles with the surrounding cell, so that particles trapped in- 
side the eddy core remain there. Under non-steady flows, particles can leak out of the 
vortices, be ejected through filamentation processes, or even the eddy can be destroyed 
i6i ( Basdevant and Philipovitchl . Il994t iHua and Kleinl , 1 199a ) . 



The OW parameter has some well-known limitations. On one side, it assumes that velocity 
gradients are slowly evolving in time, which is only valid inside relatively coherent regions. 
On the other side, this parameter fails to properly identify r egions with different mixing 



properties when eddies are stationary and have axial symmetry ( Lapevre et all 1 19991 ) . There 



have been some attempts to solve so me of these lim i tations by extending it to consider the 
time evolution of velocity gradients ( Hua and Kleinl . 1 19981 : IHua et all Il998t lLapevre et al 
19991) . 



169 Another approach to the characterization of flow structures is to assume a Lagrangian 

170 viewpoint, t hat is, to look explicitly at transport properties from the analysis of particle 



i7i trajectories ( Pierrehumbertl 1991 : Ridderinkhof and Zimmerman . 1992t Miller et al. . 19971 



Haller and YuaiJ.l200dlCoulliette and Wiggins!. l200ll :lKoh and LegrasL l2002Hd'Ovidio et all 
20041 IShadden et all 120051 iMancho et all 120061 ). In contrast to the OW method and in 
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general to Eulerian diagnostics, this approach requires knowledge of the time variability of 
the velocity field, as well as the use of an integrator for generating the trajectories. Sev- 
eral methods exist, one of the simplest being the Finite-Size Lyapunov Exponent (FSLE). 
The FSLE is a generalization of the concept of Lyapunov exponent to finite separations. 
The standard definition of Lyapunov exponent refers to the exponential rate of divergence, 
averaged over infinite time, of infinitesimally closed initial points. The FSLE technique 
keeps the original idea of capturing the rate of divergence between trajectories, but over- 
comes the limit operations. Thus, it is (and has been shown to be) rather appropriate 



to manage r e al dat a. The FSLEs were introduced for turbulent flows (jAurell et al 



1997; 



Art ale et all 119971 ) aiming at studying non-asymptotic dispersion processes. Since then, 



they have been used f or two complementary goals: for characterizing dispersion processes 



(jLacorata et all 120011 ). an d for detecting and y i sualizing Lagrangian structures (e.g. trans 



port barriers and fronts') (IBoffetta et all l200ll ; iKoli and Legrasl 120021 ; Ijoseph and Legrasl 



2002[ Id'Ovidio et all 12004). In the framework of this paper, in order to compare with the 



OW parameter we will focus mainly on the second use. 
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Several methods allow to calculate the FSLEs. In the simplest scheme, for each instant t 
and each point x, one follows in time the evolution of a tracer started in x and of another 
probing tracer located at a distance So from it. The integration is stopped when the two 
tracers have reached a final separation Sf > Sq. From the time interval, r, to reach the final 
separation, the FSLE is defined in the following way: 



\(x.,t,6 ,6f) 



1, S f 
-log/. 
t do 



(2) 



In order to reduce the dependence on the direction of the probing tracer, the algorithm is run 
choosing three points forming an equilateral triangle around x. We stopped the integration 
when any of these three points reaches a separation Sf from the trajectory started in x. 
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Maxima (ridges) of Lyapunov values are typically organized in convoluted, lobular lines. For 
the case of the backward calculations, these lines can be interpreted as the fronts of passively 
advected tracers. This interpretation can be understood in a very qualitative but effective 
way, considering that fronts typically separate fluid patches of different origins. Therefore, 
the separation in the past of a couple of points is largest when the couple is initialised 
exactly over the front. This argument can be rephrased in a slightly more rigorous way in 
the context of dynamical systems, interpreting the line-shaped regions of fastest separa- 
tion -cither backward or forward in ti me- respectively as the unstable and stable manifolds 



2001 



of the hyperbolic points in the flow (IHaller and YuanL 120001; IB offetta et all 12001b lH allcr 



Koh and Legrasl 120021 ; I Joseph and Legrasl l2002l ; fd'Ovidio et all 12004 flVIancho e~ 



;< 



20061 ). The effect on advection of hyperbolic structures is sketched in Fig.Q] Due to the con- 
vergent field along the stable manifold and the divergent field along the unstable manifold, 
a passively advected tracer is deformed as in Fig. [TJ developing a front along the unstable 
manifold and a gradient orthogonally to it. Due to the hyperbolic structure of which the 
manifold is a part, the tracer front approaches the manifold exponentially fast. The dis- 
tance S r between the tracer front and the manifold depends on the initial front-to-manifold 
distance Si, the exponent A of the manifold, and the time of integration t: 



Si exp(— At). 



(3) 



Manifolds characterized by higher exponents have therefore a stronger effect on tracers, 
shaping a front in shorter times and being more visible in tracer distributions. 
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Note that for a time-dependent velocity field, the sketch of Fig.[T]holds only if the hyperbolic 
structures evolve in time on a time scale slower than the trace r advection, so that the tracer 
can actually relax over the manifold. See lLehahn et al.l 1 20071 ) for more details on the FSLE 
computation from altimetric data. 



222 Mixing properties can also be diagnosed by Lyapunov exponent calculations, either con- 

223 sidcring the exponential separation in the fut ure as a measure of t racer dispersion, or by 

224 combining forward and backward information (jd'Ovidio et al.l . 120041 ) . Since in this work we 



aim at a direct comparison with advected tracers (SST), we will focus on the backward 
calculation and compare the location of the manifolds detected by the FSLEs with tracer 
fronts. 



228 As it is clear from Eq.([2]), the FSLEs depend c ritically on the c hoice of two length scales: 

229 the initial separation <5 and the final one Sf. Id'Ovidio et al.l (J200J) argued that Sq has 

230 to be close to the intcrgrid spacing Ax among the points x on which the FSLEs will be 

231 computed. In fact, a <5o larger than intergrid spacing would allow to sample manifolds of 

232 strong Lyapunov exponents in more than one grid points, while <5o smaller than the intcrgrid 

233 spacing w ould not allow to fol lo w a manifold on the sampling grid as a continuous line. 

234 Following Fd'Ovidio et"aD (J2004I ): iLehahn et"aH (J2007t ). in this work we have set So = 0.01° 

235 (approx. 1 km) in order to match the resolution of SST images. We have set Sf = 1° i.e., 

236 separations of about 110 Km, that is the order of magnitude of the eddies' radii detected 

237 by altimetry. Values of Sf smaller or larger up to 50% do not change significatively the 

238 calculation. The time of integration for finite-size Lyapunov exponents varies from one 

239 point of another, being small for strong values and vice versa. Inverting Eqj^l 



1, S f 

T 1o S'a ■ 
A d 



(4) 



We found typical values of finite-size Lyapunov exponents in the range 0.1 — 0.2 days 
corresponding to integration times of resp. 46 and 23 days. 



243 3 Data 



244 3.1 Sea-Surface Height 
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In the Mediterranean sea, despite the weak signal intensity and the coarse space and time 
resolution of the altimetric tracks, several studies have shown the relia bility of the alti- 
metric data to analyze its dynamics, parti c ularly in the Alge rian basin ( Vignudellil . Il997t 
Ifjouzinac et al.l . Il998t lLarnicol et all l2002t I Font et al.l . 120041) . In this study we have used 



Delayed Time Maps of Absolute Dynamic Heights (DT-MADT) produced by Collecte Local- 
isation Satellites (CLS) in Toulouse (France) specifically for the Mediterranean sea, which 
combine the signals of ERS-ENVISAT and TOPEX/Poseidon-JASON altimeters. 



Altimetric data are processed including usual corrections (sea-state bias, tides, inverse 
barometer, etc.) and improved orbits. From several corrected sea-surface height files, a 
conventional repeat-track analysis is performed to extract the Sea Level Anomaly (SLA) 
relative to a mean profile: data arc re-sampled along the mean profile using cubic splines 
and differences relative to the mean profile arc calculated. SLA along-track data arc then 



257 filtered and subsampled. The filters used are a non linear median over 3 points (roughly 21 

258 km) followed by a low pass along track linear Lanczos fi lter (with a cut-off wavelength of 42 



259 km). SLA data are then subsampled every other point (jSSALTO/DUACS User Handbook , 

260 120061 ). Finally, SLA maps are built using an improved space/time o bjective analysis meth od, 
26i which takes into account long wavelength errors, on a regular grid ( Le Traon et alj . ll998 ) of 



(1/8)° x (1/8)° every week. Then, Sea-Surf ace Heights (SSH ) are finally obtained by adding 
to the SLA a Mean Dynamic Topography ( Rio et all 120071 ) . 



264 For each data set gcostrophic velocities arc estimated as usual: 

g dh ssh g dh ssh 

jRt oq> jRrcoacp oX 



266 where h sa h is the SSH, g is gravity, / the Coriolis parameter, Rt the Earth radius, </> the 

267 latitude and A the longitude. The data analyzed spans from January 1, 1994 to December 31, 

268 2004. From this data set we study representative days and we also construct climatologies. 

269 Finally, a Runge-Kutta integrator of fourth order and a time step of 6 hours has been 

270 used to obtain backward and forward trajectories in the velocity field for the calculation of 

271 FSLEs and for the advection of a synthetic tracer. The gcostrophic velocity field has been 

272 resampled in space and time with a multilinear interpolator. 



273 3.2 Sea-surface temperature 



274 We used sea-surface temperature (SST) data from the AVHRR sensors on board NOAA 

275 satellites, downloaded from the HRPT station at the Institut de Ciencies del Mar (CSIC) 

276 in Barcelona. For the single day analysis we have chosen images for which there were good 

277 quality sea-surface temperature images presenting a large variety of eddy structures and 

278 different intensities of the Algerian current: July 9, 2003; April 7, 2004; June 30, 2004. SST 

279 images have a resolution of 1.1 km at the nadir. 



4 Results 



4-1 Eddy representation 



282 Figure [3 shows the spatial distribution of OW, which appears as a set of vorticity-dominated 

283 (W < 0, blue) regions surrounded by strain-dominated lobular structures (W > 0, red) 

284 embedded in a background field of small values of W. On the other hand, Fig. [3] shows the 

285 spatial distribution of FSLEs which has very different patterns characterized by a tangle 

286 of lines. These are the locations were FSLE are large, approximating unstable manifolds 

287 and corresponding to transport barriers embedded in a background field with A ~ 0. The 

288 unstable manifolds arise from fluid stretching and they are therefore Lagrangian analogs of 

289 the deformation OW regions. 

290 A first important difference in daily maps of OW and FSLE values is the presence of eddies 

291 in the first and of filaments in the second. In particular, in the FSLE map there are no 

292 enclosed regions and in some cases (e.g. for the eddies over the Balearic abyssal plain), the 



293 same manifold connects multiple lobes, spanning a region of several eddies. When comparing 

294 these patterns with the stream-function of the flow, SSH in this case, it can be observed 

295 that in general eddies localized by the OW parameter and the centers of the spiralling FSLE 

296 lobes are in good agreement with the extrema of the altimctric field, as one would expect. 

297 Across the Algerian current, the OW field identifies some possibly spurious eddies that do 

298 not appear in cither the FSLE or SSH. This is due to the fact that both SSH and FSLE 

299 are not invariant under a transformation of coordinates to a frame of reference moving at 

300 a constant velocity with respect to the original (Galilean transformation) and therefore, 

301 eddies are hidden or partially hidden by the presence of the Algerian current. Since OW is 

302 Galilean invariant it is able to detect these eddies but it fails in the detection of the Algerian 

303 current which appears as a coherent structure characterized by manifolds (barriers) parallel 

304 to the altimetric streamlines in the FSLE picture. 



305 Another key difference between both fields is linked to the time-evolution of coherent struc- 

306 tures. From Figs. [2] and [3] it can be seen that eddies displaying a similar size and intensity 

307 in the OW map may show very different features in the FSLE map. This is the case for 

308 instance for the western Alboran eddy and some of the mesoscale eddies in the Algerian 

309 basin (especially during July 30, 2004). For the Alboran vortex, the lobe spiral is very 

310 tight, almost resembling the concentric altimetry lines. For the eddies in the Algerian basin, 
3ii the lobes are loose, much less localized, and di scordant with the altimetric contours. It is 
312 interesting to notice that previous studies (e.g. Ilsern-Fontanet et all 120061 . and references 



313 therein) have shown that vortices in the Algerian basin propagate at velocities of the order 

314 of 5 km day -1 in contrast to the western Alboran eddy, which is almost stationary. In 

315 the ideal case of a time-independent velocity field, particle trajectories, as well as the set 
3i6 of lines in which Lyapunov exponents are organized (which approximate material lines), 
317 coincide with altimctric isolines. For the case of an eddy, they appear as concentric, closed 
3i8 manifolds. In a time-dependent flow, the identity between trajectories and altimetric iso- 

319 lines is lost and the differences quantify the time variability. The material lines resemble 

320 concentric circles for slowly evolving persistent structures like the Alboran gyres and as- 

321 sume a complex shape in the case of more dynamically active eddies, like the eddies along 

322 the Algerian current and some of the eddies over the Algerian basin. This phenomenon 

323 has an important effect on the transport properties. Due to the fact that material lines act 

324 as transport barriers, eddies with a low time variability and concentric unstable manifolds 

325 have a smal ler water mass excha nge with the surrounding compared to dynamical active 

326 vortices (see lLehahn et al.l (|20071 ) for a discussion of the eddy time variability in connection 



327 with phytoplankton pattern formation). 



328 A third important difference between the OW and FSLE approaches is the spatial scale of 

329 the structures detected. The OW parameter is bounded by the resolution of the altimetric 

330 data (1/8°) and the low-pass filters applied during the construction of altimetric maps. 

331 This limitation does not hold for the FSLE that is based on trajectory calculations, whose 

332 length scale is a combination of both space and time variability of the velocity field. Lobes 

333 and filaments below the altimetric resolution appear, especially in the more dynamically 

334 active regions, like the fast evolving eddies formed downstream of the Algerian current or 

335 the lobes south of the Balearic islands. Very thin filaments are also associated with the 

336 Algerian current. 



337 Jf.,2 Detection of tracer fronts 



338 lf.,2.1 Synthetic tracer 

339 First, we test the ability of OW and FSLE to characterize tracer distribution in an ideal but 

340 realistic situation: three sets of particles are advected by the velocity field estimated from 

341 altimetry. To this end we put three sets of particles distributed on a square grid centered 

342 over three different dynamical structures: the slowly evolving Alboran eddy, the Algerian 

343 current, and the strongly interacting eddies in the easternmost part of the Algerian basin 

344 (dashed boxes in Fig. HJ. Particles in the eddy regions have been advected for two weeks. 

345 A shorter advection time of one week has been used for the particles initially placed over 

346 the Algerian current, due to the strongest velocity field in this region. The advection time 

347 that we have chosen is such that the particles are kept close to the dynamical structures we 

348 want to study: a larger advection time does not change the results discussed, but increases 

349 the dispersal of the particles over several structures. 

350 The tracer released over the west Alboran eddy shows a regular, circular pattern well corre- 

351 lated with both the OW and the FSLE maps. This pattern differs from the more deformed 

352 distribution for the tracer released over the Algerian basin. In this region, the tracer ap- 

353 pears spread on several eddies connected through thin filamentary structures. The bound- 

354 aries detected by the OW parameter (we plot in Fig. [4] top, the lines W = 0) provide 

355 an approximate picture, often underestimating the size of the eddy cores and providing no 

356 information of the patterns followed by the tracer exchanged from one eddy to another. 

357 In contrast, the FSLE lines of intense stretching reproduce with remarkable accuracy the 

358 tracer eddy boundaries, as well as the tangle of spiralling filaments that connect them. This 

359 is seen in Fig. [4] (bottom) where we plot the regions where the FSLE has values larger than 

360 0.2 day^ 1 . As discussed above, these regions are essentially one-dimensional lines. They 
36i behave as material unstable manifolds of the advecting flow: they are almost perfectly lo- 

362 cated along the tracer boundary. Small deviations may be attributed to different reasons: 

363 (i) a tracer front approaches the manifold exponentially fast, so that a residual distance 

364 remains for a finite time of integration; (ii) only the manifolds with an intensity larger than 

365 0.2 day -1 have been plotted and other weaker lines may also act as transient transport 

366 barriers. The residual distance from a manifold can be estimated by Eq. [3J considering an 

367 initial tracer to manifold distance of about 100 km, an exponent of 0.2 — 0.3 day^ 1 for 

368 the manifold (Fig. [3]), and an advection time of 15 days, we get a separation between the 

369 manifold and the tracer front of a few km. Note that the thin filament at 7°W, 37.5°N is 

370 below the altimetric resolution, and appears in both the tracer advection and the FSLE 

371 map. 

372 Thin filaments also appear for the case of the tracer released over the Algerian current. 

373 The OW is designed to detect vortices and therefore its use for barrier detection along the 

374 Algerian current is, strictly speaking, improper. Nevertheless, the gyre due to the mesoscale 

375 eddy located at 0.5°W, 36.5°N is correctly predicted. Not surprisingly, features due to 

376 smaller and rapidly evolving eddies (like the ones located along the current) , as well as the 

377 barrier effect due to the current itself, are not detected. Interestingly, an OW signal appears 

378 in correspondence to the Almeria-Oran front (well visible in the southern boundary of the 

379 tracer), probably as a signature of the secondary Alboran eddy. Such a signature is composed 

380 of broken and irregular structures, but nevertheless is in phase with the southern front of 
38i the tracer. The OW does not provide any indication of the filament intruding the Algerian 

382 basin. In contrast with the lobular structures detected for eddies, the FSLE shows for the 

383 Algerian current meandering lines that follow the African coast. Such lines are in almost 

384 perfect agreement with the tracer distribution, well indicating the region of intrusion in 
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385 the Algerian basin. A spiralling lobe at the location of the eddy also detected by the OW 

386 method at 0.5°W, 36.5°N is correctly localized. The Almeria-Oran front is well detected 

387 by a manifold that follows uninterrupted the Algerian current, marking a transport barrier 

388 that confines the water masses coming from the Alboran sea to an isolated tongue along the 

389 African coast. In agreement with this picture, the tracer has no intrusion with such coastal 

390 water, being initialized northern to such a manifold. 



391 J^.2.2 An observed tracer: SST 

392 Figure O depicts the temperature distribution for the three days selected and Fig. [5] shows 

393 some zooms corresponding to the dashed squares in Figs. [51 [3] and O 

394 Figure [B^ shows the situation of a relatively isolated westward propagating Algerian vortex. 

395 The center of the eddy is well located by the OW parameter, but the strongest SST gradients 

396 are beyond the outside of the W = which identify the vortex core. FSLE lines reproduce 

397 the two-lobe structure of the SST positive anomaly contered in IE, 37. 5N, although the 

398 size is overestimated. Furthermore, south of this vortex there is a filament of colder waters 

399 that approximately follows the transport barriers depicted by FSLE. However, this example 

400 also shows one of the limitations of the approach: FSLE and OW strongly depend on the 

401 quality of altimetric maps. A coastal eddy in the Algerian coast is clearly observable in the 

402 SST images but is not properly captured by SSH, and therefore partially missed by OW 

403 and FSLE. 

404 Figurc[Bb depicts the Alboran sea in spring (April 7, 2004). In its westernmost part, close to 

405 the Gibraltar strait, a water mass of warm waters, surrounded by colder waters, is trapped 

406 within the western Alboran gyre. Both the OW and the FSLE maps show a regular and 

407 circular barrier although both seem to underestimate the radius of the object, maybe due 

408 to the location of the vortex in altimetric data. The Almeria-Oran front on the eastern part 

409 of the image is observable in both SST and FSLE but not in the OW field. Proceeding to 

410 west a large eddy attached to the coast is quite properly identified in all fields: SST, FSLE, 
4ii SSH and OW. However, in the middle of the image, where the eastern Alboran eddy is 

412 usually located there is a poor coincidence between patterns calculated from altimetry and 

413 the SST image. In particular, the SST fronts observed at 3°W are almost perpendicular to 

414 the FSLE lines there. 

415 Figurcint shows an example of the strong signature of the Algerian current. In contrast with 
4i6 the previous case, the FSLE map reproduces with great accuracy the SST distribution. The 
417 dynamical barriers due to the presence of the jet along the coast are parallel to the SST 
4i8 isolines. In analogy to what was observed for the synthetic tracer, the jet has no signature 

419 on the OW parameter. On the other hand, several mesoscale eddy boundaries with no effect 

420 on the SST pattern also appear in the OW image. Centered at 37°N 2°E there is a large 

421 coastal eddy probably generated by the destabilization of the Algerian current. This eddy 

422 is clearly identified in all fields. However, the most remarkable pattern is the deflection 

423 of coastal waters from coast to the open sea due to the presence of this eddy, which is 

424 observable in FSLE as well as in SST. Note also that water masses in the inner part of the 

425 vortex are bounded by transport barriers and therefore trapped within the eddy. 

426 Finally, Fig. [BJi shows the easternmost part of the Algerian basin where several eddies 

427 strongly interact. Notice that this is one of the situations analyzed using the advection of 

428 ideal particles in the previous section. As expected, the match between real SST data and 

429 the FSLE is as good as before. The manifold tangle observed in the FSLE map can also 

430 be seen in the SST image and the cores of the eddies on which these manifolds wind are 
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431 identified also by the OW parameter. An isolated region along the Algerian coast (7E, 37N), 

432 appearing in the SST image as a strong cold anomaly, is also fairly well identified by FSLE 

433 and not by the OW parameter. 



434 4-3 Time averaged fields 



435 In order to generalize the comparison performed on instantaneous cases, we focus now on 

436 the relationships between spatial field distributions averaged over time. To this end OW 

437 and FSLE have been computed for 10 years of data (1994-2004). The OW parameter is 

438 calculated for each altimetric image (at one week time resolution) while an FSLE map is 

439 generated each two days. 



440 4-3-1 Distribution of barrier-type lines 

441 First, we estimate the fraction of time during which each spatial point is visited by barrier- 

442 type lines. By barrier-type lines we mean lines (ridges made of local maxima in FSLE and 

443 the W = isolines in OW) which could be interpreted, at least during some short time, 

444 as a transport barrier. To this end we count how many times each grid point pertains to 

445 one of such lines and divide by the total number of observations for that point (the length 

446 of the time-series). Figure [7] shows these local probabilities of having a transport barrier 

447 estimated using OW and FSLE (and denoted by Pw and P\, respectively). In agreement to 

448 what we found in the previous sections, we observe that the pattern corresponding to OW 

449 parameter has a more regular, patchy structure due to the presence of Eulerian eddies, and 

450 no clear signature for features that are not dependent on individual vortices, like the Algerian 

451 current which appears instead in the FSLE map. This suggests the idea that features that 

452 are common in both maps might be associated to long-lived non-propagating vortices. The 

453 most evident example is the western Alboran eddy. Other examples are the pattern observed 

454 around longitude 4°E and 38°N whi ch is associated to the region of vortex detachment from 

455 coast discussed in previous papers ( Isern-Fontanet et all 120061 ) . the pattern observed east 



456 of Eivissa island (1.5°E and 38.5°N) or the Almeria-Oran front. 

457 After observing that OW and FSLE provides very different estimations of the propensity of 

458 the points to belong to a barrier-like line, the next question is to identify which one is closer 

459 to give a true barrier intensity and location. As outlined in the introduction, a characteristic 

460 of the incoming Atlantic waters with respect to the resident Mediterranean ones is their 
46i lower salinity. If there are significant barriers to the spread of this surface water into the 

462 Mediterranean, climatological distributions of salinity should locate them. Figure [8] shows 

463 the climatological salinity obtained from MEDATLAS-II data set. When comparing Figs. 

464 [7] and [8] it is evident that the estimation computed from FSLE provides a better picture, 

465 although the smoothness of the climatological salinity makes difficult the comparison. On 

466 Alboran, salinity increases eastwards. Close to the Almeria-Oran front the isohalines are 

467 almost aligned width the front, as depicted by FSLE. Eastwards, the intrusion of fresh 

468 waters along the coast following the Algerian current and the associated intense northwards 

469 gradient matches quite well the patterns observed in Fig. [Jj (bottom). Proceeding to the 

470 cast, at the entrance of the Sardinia channel FSLE barrier density is concentrated close to 

471 the coast, in correspondence with low salinity waters being also confined to this area. In 

472 the middle of the basin, the probability of having transport barriers is more homogeneous 

473 due to the propagation of Algerian eddies and consequently climatological salinity is more 

474 uniform. 
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475 4-3-2 Time averages of FSLE and other Eulerian diagnostics 

476 Results in previous sections suggest that high FSLEs are preferentially located on ener- 

477 getic structures like at the eddy periphery or along jets. This suggest a possible correlation 

478 between the FSLEs and other Eulerian diagnostics different from OW, such as the eddy 

479 kinetic energy (EKE) or the strain. In fact, in a recent study of the Tasman sea, it has 

480 been shown that averages of both the EKE an d the strain have a compact relationship 
48i with averages of finite-time Lyapunov exponents ( Waugh et all 120061 ). Although such rela- 



497 



504 



482 tionships are valid only on time average, and therefore cannot be used for the detection of 

483 instantaneous, local transport barriers, the correlation between EKE, strain, and Lyapunov 

484 exponents w ould provide a simpl e way to estimate Lyapunov exponent climatologies, as also 

485 suggested in I Waugh et alj (2006). Finite-time and finite-size Lyapunov exponents only differ 



486 on the determination of the integration time (set a priori for finite-time, given indirectly by 

487 prescribing the final separation of the trajectories for finite-size). For this reason, we tested 

488 the hypothesis above by computing for our data period time averages of strain rate and of 

489 EKE. 

490 Climatologies of FSLE, EKE, and strain rate are shown in Fig. [9] Similar structures appear 

491 in regions containing persistent objects, such as the Alboran gyre and the Algerian current. 

492 A large-scale gradient is also observed, with more signal in the south than in the north of 

493 the Algerian basin, and positive anomalies in the Sardinia channel. However, there are more 

494 localized mesoscale anomalies in the FSLE field. 

495 In Fig. [TO! joint distributions of FSLEs vs strain rate and EKE are presented. Although 

496 theres is a general positive correlati on among these qua ntities, they show a much looser 
relationship than the one observed by ( Waugh et al.l . l2006l ) at lower resolution and on shorter 



498 integration times. The lack of compact relationships on the scatter plots, as well as different 

499 distribution of anomalies among the various diagnostics is a clear confirmation that the 

500 Lyapunov exponent provide a complementary information in respect to Eulerian diagnostics 

501 even in a climatological sense, when the Lyapunov calculation is done at high resolution. The 

502 comparison with other Eulerian diagnostics (like the vorticity and the v orticity gradient , 



503 that is sometimes also used for transport barrier detection, see for instance (jPaparella et al 
19971 )) confirm these observations. 



505 5 Summary and discussion 



506 In this work we have exploited the rich mesoscale activity of the Algerian basin for comparing 

507 Eulerian and Lagrangian diagnostics of transport under various conditions. This comparison 

508 unveiled two coexisting pictures: an Eulerian view, dominated by mesoscale eddies, and 

509 a Lagrangian view, characterized by interconnected lobular structures and submesoscalc 

510 filaments. The first picture describes the eddies that populate the mesoscale turbulence; 
5ii the second one describes the spatial structures of tracers advected by these eddies. Eddy 

512 cores are usually identified by setting a threshold (W = 0) on the OW parameter, and this 

513 W = line may be a rough candidate to a 'transport barrier' since eddy contents are seen 

514 to remain coherent in several situations. More precise candidates to transport barriers are 

515 obtained from local maxima (ridges) of the FSLE computed backward in time. 

5i6 As diagnostics of transport barriers, both techniques are based on heuristic arguments 

517 related to the time variability of the velocity field. The OW use is based on the consideration 

5i8 that for time-invariant fields, Lagrangian and Eulerian barriers coincide. For this reason, 

519 with the OW method no structure typical of time varying fields like filaments and lobular 



13 



520 patterns can be detected but nevertheless remarkably good results can be obtained for 

521 stationary eddies, in our case the first Alboran gyre and some of the most persistent eddies 

522 over the Balearic abyssal plain. The FSLE calculation goes one step further, assuming quasi 

523 stationarity not for the velocity field, but for the lines of maximal divergence (unstable 

524 manifolds of hyperbolic points). When this time-scale separation between the Lagrangian 

525 structures and the advectcd tracer holds, the tracer front relax exponentially fast over 

526 the ridges of large Lyapunov exponents that therefore mark the front position. Another 

527 advantage of the FSLE method is that the effect of both eddies and jets are remarkably 

528 dealt with. 

529 We have tested the barrier candidates against a synthetic tracer advected by the geostrophic 

530 field. This test sets an upper limit on the confidence on FSLE and OW assumptions, since it 

531 excludes the effect of any factor acting on the tracer except the geostrophic velocity field. At 

532 the meso- or larger scales, the OW was found to locate correctly (with errors of tens of km) 

533 tracer patches when they were well confined inside eddies. Fluid parcels intruding eddies's 

534 along thin filaments arc not captured however. Such filaments, as well as the submesoscalc 

535 structure of tracer fronts, were found from the FSLE calculations. 

536 In fact, the advection of a synthetic tracer can be considered itself as a direct, Lagrangian 

537 diagnostic of transport barriers. This approach however has some limitations. Lagrangian 

538 structures are sampled unevenly, since the tracer tends to spend longer times in regions with 

539 low velocities. For this reason, even strong barriers may require a carefully choice of the 

540 initial conditions and of the integration time. An example of this appears in Fig. 2J where 

541 part of the Almerian-Oran front is not entirely visible and would have required to probe 

542 the velocity field with a larger number of tracer blobs. When doing the tracer experiment of 

543 Fig. [J] we also observed that, not surprisingly, the tracer has a long residence time over the 

544 Balearic Abyssal Plain and therefore for long integration times most of the time the tracer 

545 shades the barrier of this region only, independently of its initial condition. On the other 

546 hand, a reduction of the integration time increases the gap between the tracer front and 

547 the actual barrier position, since the tracer front has a shorter time to relax to it (see the 

548 discussion of Eq. [3|. In fact, an attempt to overcome these limitations requires the use of 

549 backward trajectories (for which barriers are regions of maximal separation) and of variable 

550 integration times, basically yielding a FSLE-type calculation. For this reason, the FSLE can 

551 be considered an optimised way of advecting a passive tracer for the detection of transport 

552 structures. Note also that, besides guaranteeing a uniform sampling of the velocity field and 

553 the precise localization of transport barriers, the FSLE also provides the information of the 

554 intensity of the barrier at the same cost of tracer advection. 

555 The representation of advection in the Algerian basin however is based on altimetry and 

556 therefore is limited by the lack of representation of any non-gcostrophic component of the 

557 velocity fie ld as well as unreso lved or poorly resolved structures due to space and time 

558 resolution ( Pascual et all [2006J), and noise. The comparison with SST images allows to 



559 verify that such effects are not strong enough to spoil the FSLE usefulness in locating 

560 Lagrangian features. Note that the filament tangle connecting different eddies as well as the 
56i spiral structures inside eddy cores cannot appear in any way in a time-independent field 

562 (i.e., by using a snapshot of the velocity field), where one can show that orbits correspond 

563 to altimetric contours, and eddy barriers to closed contours. Nevertheless, the slow time 

564 variability of a few altimetric frames is enough to generate Lagrangian structures which 

565 correspond well to SST images. We also expected an effect of altimetric errors especially 

566 for the structures close to the coast, but with the exception of the position of the first 

567 gyre in the Alboran sea and in spite of the chaotic nature of the geostrophic advection, the 

568 localization of features from Lyapunov exponent maps appears quite reliable. 

569 Climatologies of several Eulcrian diagnostics and cither barrier densities or average values 



11 



578 
579 
580 
581 
582 
583 
584 



of Lyapunov exponents also showed distinct different patterns. Given the fact that the main 
differences between Lyapunov and Eulcrian diagnostics are related to small scale filaments 
and in the temporal variability of mcsoscale turbulence, it is not difficult to expect different 
relation between Eulerian diagnostics and Lyapunov exponents computed at much larger 
scale and on shorte r time periods. This is the case of the 0.5 deg. and 15 days calculation 
of in IWaugh et al.l ( 200(a ) (resp. 50 times larger and 3 times shorter than in our case). 
We illustrate this in Fig. 111! where the Lyapunov exponents have been computed with 
initial separations of 0.5 deg. for the same day of Fig. [3Jt. Indeed, at the resolution of 
the Mediterranean altimetry product (1/8 deg.) not even the patterns of two Eulcrian 
diagnostics like the strain and the eddy kinetic energy appear to be correlated (Fig. [5] b 
and c) . An interesting exercise would be the calculation of correlation between eddy kinetic 
energy and Lyapunov exponents computed at various scales, in order to see whether there 
exists a critical lengthscale or a smooth transition. In fact, the finite-size Lyapunov exponent 
has been introduced in turb ulence theory with the aim of probing a velocity field at various 



scales. General arguments (jAurell et al.l . 119971 ) suggest the presence of a cross-over scale 



in the Lyapunov calculation separating an exponential regime from other diffusive regimes 
at the scale of the smallest resolved eddy. For the case of altimetry, this cross-over can be 
expected to appear at spatial scales of the order of 1 deg. 



588 
589 
590 
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Note that the Lyapunov exponent is not an advanced replacement of Eulerian diagnostics, 
but a complementary tool incorporating tracer information. For instance, the association of 
lobes with eddies that in some cases we observed only holds for stationary eddies. If used 
systematically as a Lagrangian eddy detection, this approach would lead to detect spurious, 
small-scale spatial signal, that in fact comes from the temporal variability of large-scale 
eddies. The OW and FSLE method are also complementary from a technical viewpoint. The 
OW method is a differential technique. It is thus affected by small-scale noise and especially 
by discontinuities. These kind of features indeed can appear in the presence of eddies with a 
size close to the spatial resolution. On the contrary, the FSLEs are computed by integrating 
particle trajectories and thus are not strongly affected by local noise or discontinuities. The 
complementarity of the information detected by FSLE in respect to the OW parameter and 
other Eulerian diagnostics is also emphasised in the climatologies of Fig. [7| that showed how 
the differences detected in single day analyses are indeed a persistent feature. 
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The main drawbacks of the FSLEs are a relatively increased complexity in the calculation 
(with respect to the OW) and especially the need of extended velocity data sets, both in 
time and in space. This can be a problem if the observational region is limited (either in 
time or space) , since some trajectories may go out from the bo undaries (this is in fac t at 
the basis of other methods used to characterize flow structures ( Schneider et all 120051 )). 
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The OW on the contrary can be calculated from a snapshot of the velocity field. Indeed , 
its easiness allows t o compute it even from in situ d ata (e.g. IStocker and Imbergerl . 120031 ; 
Testor and Gascard 120051 ; iPallas and Viudea . 120051 ) and then obtain a first approach to 



the flow topology and its effect on the transport. Furthermore, although OW is a dif- 
ferential technique and therefore more sensible to velocity errors or discontinuities, it is 
possibl e to approach i t using kinetic energy at least for the identification of vortex core 



edges ( Paparella et all 119971 : llsern-Fontanet et all 120041 ). If the vortex is assumed to be 



close to axial symmetry then the OW formula can be written as 



1 fdEOr) 

614 W ex - —^ J - 

r V or 



(6) 



where r is the radius from the vortex center and E(r) its density of kinetic energy at that 
distance from the center. This implies that the line of W — will correspond to a local 
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617 maximum of energy. 

6i8 The comparison with SST and synthetic tracer patterns also shows that the errors in the 

6i9 OW method and the FSLE technique are of very different nature. The OW fails with both 

620 the real and synthetic tracer mainly for detection of the filamcntation process, due to its 

621 intrinsic lack of representation of the time variability of the velocity field. In contrast, the 

622 FSLE is limited not by its assumption of time-scale separation between barrier and tracer, 

623 but because of the limits of the altimctric data. This is clearly shown by the fact that the 

624 FSLE performs almost ideally for the synthetic case while only providing a more qualitative 

625 estimations of SST fronts. In perspective, envisaging the availability of assimilation and 

626 observational velocities field of increased quality and spatiotcmporal resolution, one can 

627 expect a decrease in effectiveness of the OW technique in favour of the FSLE (or any other 

628 Lagrangian technique) at least for the detection of (sub-)mesoscale transport barriers. 



629 6 Conclusions 



630 The advantages of a Lagrangian tool when compared to Eulerian diagnostics in capturing 

63i tracer structures would be trivial for model data. What is not trivial is that altimetry 

632 data are amenable to such a comparison, given that the differences arise in presence of 

633 tracer filaments, i.e., at a scale similar or even smaller than the nominal resolution of 

634 altimetry data (and much smaller than the cutoff in the altimetry energy spectrum). At 

635 this scale ageostrophic components (unresolved by altimetry data) as well as SST active 

636 dynamics is supposed to play an important role. The fact that ncvcrthclss altimetry-derived 

637 patterns approach SST structures means that these unresolved compone nts of the flow 



638 are e ither of weak intensity or act in phase with the horizontal stirring (jLehahn et al 

639 



20071 ) . The recent availabili ty of realistic, submesoscale resolving biogeochemical models 



640 (e.g. (JResplandv et all 120081 )) should give the possibility of exploring the role of unresolved 



64i ageostrophic components and of tracer activity on filament detection. 

642 The fact that altimetry data can be pushed at subgrid scales when re-analised by Lagrangian 

643 tools, providing systematic information of filament location, opens several interesting per- 

644 spectives. For the case of high-resolution biogeochemical measurements, the Lyapunov ex- 

645 ponent provides an unpaired tool helping to interpret the measurements in terms of the 

646 precise position of mesoscale and sub-mesoscale fronts. The use of Lyapunov calculation 

647 can be even concieved in real-time, for targctting frontal structures during in-situ surveys. 

648 This possibility comes from the fact that the Lyapunov calculation only requires the past 

649 history of the velocity field, and from the availability of altimetry products with 1-2 weeks 

650 of delays (e.g. from AVISO). 

65i A second interesting possibility is the use of the Lyapunov exponent together with high- 

652 resolution tracer measurments for the validation of a velocity field. An improved altimetry 

653 product (as what is expected from the recent launch of Jason-2 and incoming altimetry 

654 satellites) should bring a velocity field at better spatial and temporal resolution. The pos- 

655 sibility of computing with Lyapunov calculation the position of fronts and of comparing 

656 them with high-resolution SST or chlorophyll images is one of the very few possibilities for 

657 testing the ability of the new products of better representing mesoscale and sub-mesoscale 

658 processes. This validation would be especially useful for choosing the optimal combination 

659 of spatial and temporal resolution of surface current products. A similar approach can be 

660 used also for validating assimilation products and in general for tuning the spatiotcmporal 
66i resolution of circulation models. 
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662 Given the fact that transport barriers in the ocean evolves on a much slower time scale than 

663 advected tracers, one can also approximate the short-time advection of a released tracer as 

664 a relaxation dynamics towards the stronger transport barriers and in this way estimate the 

665 tracer ditribution in the near future. The most evident application of this system would 

666 be, for example, the prediction of the e volution of pollutan ts. Other Lagrangian techniques 



have already approached this problem (JLekien et al.l . 120051 ) 
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Fig. 2. Okubo-Weiss parameter (W, coded in colors) computed from sea-surface 
height (black line) fields corresponding (from top to bottom) to July 9, 2003; April 
7, 2004; June 30, 2004. The dashed boxes represent the regions over which we 
perform a comparison with SST (see Figs. [5] and [6|). 
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Fig. 3. FSLE (A) computed from sea-surface height fields corresponding (from top 
to bottom) to July 9, 2003; April 7, 2004; June 30, 2004. The dashed boxes represent 
the regions over which we perform a comparison with SST (see Fig. [5]). 
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Fig. 4. Distribution of tracer particles released over the Algerian current, the west- 
ern Alboran eddy, and the easternmost part of the Algerian basin (dotted boxes). 
Tracers in the Algerian current are advected for one week, and the other two sets 
are advected for two weeks. The release dates are chosen to have for the three cases 
the coinciding final time (June 30, 2004) which is shown in the figure. Zero-lines of 
OW corresponding to the final day are overimposed (top) as well as the line-shaped 
regions where the backwards FSLE values are larger than 0.2 days -1 (bottom), 
approximating unstable manifolds of the chaotic flow. 
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Fig. 5. General views of the sea-surface temperature fields corresponding (from top 
to bottom) to July 9, 2003; April 7, 2004; June 30, 2004. 
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Fig. 6. Comparison of OW, FSLE and temperature distribution corresponding to 
the regions shown in figures [2][5j Left column: SST (coded in grey levels) with SSH 
(thin black line) and lines of zero OW (thick black line). Right column: SST (in 
grey levels) with regions where FSLE is greater than 0.1 day -1 (red line-shaped 
regions). From top to bottom the dates are: July 9, 2003 (a and b); April 7, 2004 
(c) and July 2, 2004 (d). 
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Fig. 7. Fraction of time during which each spatial point is visited by the lines W = 
{Pwi upper panel) and by the ridges (local maxima along some direction) of FSLE 
(P\, lower panel) for the period 1994-2004. 
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Fig. 8. Climatological distribution of salinity in the area of study at 5 m depth from 
the MEDATLAS-II data set. 
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Fig. 9. Time-averaged FSLE, strain rate and eddy kinetic energy (1994-2004). 
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Fig. 10. Joint distributions of Finite-size Lyapunov exponents vs. strain rate values 
(left), and Finite-size Lyapunov exponents vs. eddy kinetic (right) from the spatial 
distributions time-averaged over the period 1994-2004. 





-5 -4 -3 




-2 




-1 







1 




? 




3 




•1 




5 




6 




7 


3 




9 


































1 


































1 

CO 1 
































































ro 




































l-~ 




























l-O 


















P 
















































1 


ro 

































































-5 -4 -3 -2 -1 



2 3 4 5 




|0.6 

0.5 

0.4 

= 0.3 

= 0.2 

0.1 

0.0 



Fig. 11. Calculation of the FSLE for July 9, 2003 with initial separation 5q = 0.5°. 
Compare with Fig. [3] a. At this resolution, Lyapunov exponents cannot be linked to 
tracer filaments. 
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